Quantum ferroelectrics near weak first order transition 
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Considering a recent experimental indication that certain doped dielectric materials, such as 
SrTi03 and KTaC>3, are on the verge of a first order quantum phase transition, the effects of order 
parameter fluctuations in those materials near such a transition are discussed. We start with a 
semi-microscopic action, where transverse paraelectric(or optical mode) fluctuations are controlled 
by hydrostatic pressure via strain coupling, which above a certain value of strength, kills the second 
order transition. Using self-consistent parquet re-summations, which retain the leading fluctua- 
tion effects in four point vertex function, we write the free energy of the system. This fluctuation 
re-normalized free energy is then used to discuss the possibility of first order transition at zero tem- 
perature as well as at finite temperature. An asymptotic analysis of fluctuation corrected vertices, 
with strain induced long range interaction, predicts small but a finite jump in the order parameter 
near the mean field quantum critical point at zero temperature. In the case of finite temperature 
transition, near quantum critical point such a possibility is found to be extremely weak. The for- 
malism is quite general and is applicable to certain class of systems undergoing weak quantum first 
order phase transition. 

PACS numbers: 64.70. Tg, 63.70.+h, 77.80. Bh 



I. INTRODUCTION 



The static dielectric susceptibility in the materials like 
SrTiC>3 and KTa03 docs not show any divergence down 
to OK, it rather saturates at a very high value (O(10 4 )). 
Ferroelectric transitions in materials which are similar 
in structure, for example in BaTiOs, are of displacive 
type and are well described by softening of a zone cen- 
ter transverse optical mode at some finite temperature. 
This scenario, that is, the absence of transition and the 
saturation at high dielectric susceptibility, has been at- 
tributed to the smallness of the zero temperature gap 
in the corresponding zone centered optical mode. Thus 
quantum fluctuations become relevant in these materials 
at low temperature. These materials are called quan- 
tum paraelectricsi. There is a revival of interest, partic- 
ularly regarding a quantum phase transition, effects of 
the quantum critical point at finite temperature^ and 
the nature of phase transition in these materials^. The 
effects of pressure^ and impurities on the dielectric sus- 
ceptibilities have been well studied experimentally. A 
general consensus arising form these experiments is that 
the application of hydrostatic pressure moves the system 
away from criticality and the possibility of phase transi- 
tion is suppressed. One needs to apply, what is termed as, 
a negative pressure to induce phase transitions in these 
materials. One way to simulate negative pressure is to 
put non-polar impurities which creates local pressure de- 
ficiencies. In this context, experimentally^ one finds that 
T c ~ (n — n c ) 2 (where n is the average impurity concen- 
tration and n c is the critical value, typically 33%) which 
matches well with the theoretical estimated^ transition 
temperature for pressure induced transition. In this case, 
the mean field T c ~ (jp + p c ) 5 (where p is hydrostatic 
pressure and p c is the critical value) . The exponent half 



is obtained when thermal fluctuations are treated at the 
Gaussian level. The similarity between effects of pressure 
and the impurity in the transition temperature can be at- 
tributed to the high density of impurity concentrations 
in these cases. Here the disorder effects seem to be small 
and a non-polar impurity essentially induces an inter- 
nal pressure. One recent experimental report^ indicates 
that these materials show phase separation near quan- 
tum critical point. This is clearly an indication of a weak 
first order transition. Unlike a strong first order transi- 
tion, order parameter fluctuations are more important 
in materials undergoing very weak first order transition. 
One needs to study the effects of order-parameter fluc- 
tuations in four-point vertices in a proper self-consistent 
scheme. In this work, we use self-consistent parquet ap- 
proximation to take proper account of the leading order 
fluctuation effects in a quantum paraelectric near such a 
weak first order transition. We consider a displacive type 
ferroelectric in its quantum limit. In these materials the 
atomic dipoles are believed to interact among themselves 
with a long range dipolar interaction and their collec- 
tive excitations give rise to various optical modes. The 
long range nature of the dipolar interaction makes lon- 
gitudinal optical modes rigid and therefore they do not 
contribute to the low temperature phenomena. On the 
other hand transverse modes do not gain extra rigidity, 
but acquire dipolar anisotropy. Thus we start with a 
model action with short range interaction without dipo- 
lar anisotropy and try to explore the results of order 
parameter fluctuations with an effective long range in- 
teractions among them, mediated by strain fluctuations. 
We also assume that the system has a strong crystal 
anisotropy, so that the two transverse directions do not 
mix up and the system splits up in two one component 
Ising type systems in two transverse directions. Since 
ferrolcctric systems do possesses strong cubic anisotropy, 
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such an assumption does not alter the basic nature of 
the transition mechanism in these systems. We there- 
fore, use an effective Ising type model and look for the 
fluctuation effects in four point vertices. To retain the 
leading fluctuation effects in the vertex function along 
with their dependencies on the non-zero polarization, we 
calculate the free energy using a parquet diagrams. This 
fluctuation renormalized free energy indicates possibility 
of a first order transition at zero temperature and it also 
gives finite temperature corrections. 



It is assumed that < tp(q, oj) >— 0. In this approximation 
the actionpQ can be written as 
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II. SUMMARY OF THE MEAN FIELD 
ANALYSIS 



The free energy is defined as the functional integral over 
Aty,P), i.e. 



Considering the system to be at a constant pres- 
sure, we integrate out the strain fluctuations completely. 
Our starting effective action describes only polarization 
fluctuations with strain induced long range interactions 
among them and takes the following form. 
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The value of P is to be determined by minimizing the 
free energy F. Stability of a thermodynamic system re- 
quires the free energy to be positive. In a Landau theory, 
which neglects fluctuations completely, the the require- 
ment becomes (u+v) > 0. In that case the system should 
undergo a second order transition. On the other hand, 
for (u + v) < 0, one needs higher order term with posi- 
tive coefficient to ensure stability. In this case the system 
should undergo a first order transition. Though a strong 
first order transition is well described by mean field the- 
ory, a theory of weak first order transition in the limit 
(u + v) — > 0~ needs proper account of the fluctuation 
effects. 



Here (j) q (uj n ) describes the Fourier transform of local 
transverse polarization </>(x,t), q is the field momentum, 
uj n = 2irn/f3 is the Matsubara frequency for Bosonic ex- 
citations, r and u are the quadratic and the anisotropic 
short range quartic coupling respectively and v is the 
isotropic long range part of the quartic coupling in- 
duced by strain and L d is the system volume in <i-spatial 
dimension. Hydrostatic pressure, as well as the non- 
polar impurity, couples to the optical mode via strain. 
It shifts the bare quadratic and quartic coupling ro to 
r = ro(l +p/po), where po is a constant. Strain fluctua- 
tion induces a long range attractive interaction between 
the dipoles and is denoted by the effective quartic cou- 
pling v. <j)q(u)) can be decomposed into two parts, P, 
the static(mean field) and ^(q,uj), the fluctuating part 
as follows, 



<t> q = P6(q,u>)+i)(q,u;). (2) 



III. FLUCTUATION CORRECTIONS TO THE 
FREE ENERGY AT ZERO TEMPERATURE 

To describe the effects of the order parameter fluctu- 
ations near a weak first order phase transition, we need 
to calculate the fluctuation renormalized four point ver- 
tex functions. Here we propose a quantum generaliza- 
tion of the work done by Gadeker and Ramakrishnan&2 
to calculate the renormalized four point vertices in par- 
quet approximation. The technique involves summation 
of leading order diagrams at each order. It has also been 
shown to be equivalent to the renormalization group^ 
calculation. It turns out that, at or near d = 4 the di- 
agrams which can be cut into two parts by cutting two 
lines in the same bubble, contribute to the leading order 
corrections. For example, the diagrams in figure PQ are 
important for the fluctuation corrections in short range 
vertices in parquet approximation. Since in this case the 
system acquires a small but non zero polarization P even 
above the mean field critical point, the polarization fluc- 
tuation becomes gaped, with the gap being proportional 
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to P 2 . Thus the free propagator for this system is given 
by, 

G' 1 (q,uj n ) = r + cq 2 + ui 2 + (u/2 + v/6)P 2 (5) 
In the paraelectric phase P = 0. A bare theory pre- 




FIG. 1: Parquet diagrams for the fluctuation corrections to 
the short interaction vertex u. It takes care of the loop cor- 
rections of the above type to all orders. Here each line corre- 
sponds to the propagator given by the equation [5], 

diets that at T = 0, the static dielectric susceptibility 
X(0,0) - G(0,0) ~ 1/r. Thus r = is a point of 
instability in the paraelectric phase, and in a situation 
where there is no discontinuity in the order parameter 
at that point, it can be identified as a quantum critical 
point(QCP). In the vicinity of the critical point, correla- 
tion length becomes large, leading to dominance of the 
order parameter fluctuations. Moreover when + w 0, 
fluctuation corrections to the four point vertices become 
important. We will try to discuss the effects of fluctua- 
tions in the vicinity of the limit |it + v\ — > 0, in devel- 
oping non-zero spontaneous polarization near the transi- 
tion point. Since in the case of a week first order tran- 
sion, initially correlation length grows significantly, the 
bare vertices get strongly renormalized. Now using the 
propagator for the polarization fluctuation (cqn.[5]) we 
find the the leading order contribution from the first two 
diagrams of the figure pQ to the renormalization vertex 
function as follows, 
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At T = 0, the frequency summation becomes an in- 
tegral and the combination of the frequency sum and 
the d-dimensional integral can be replaced by a d + 1- 
dimcnsional integral. Hence, 
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where K d is a constant which is related to the surface 
area of a c?-dimcnsional sphere of unit radius. It is given 
by, 



= (2 
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The 



= g^r ana in d=4 and d=3 respectively, 
corrections to four point vertex 5u has a logarithmic di- 
vergence as r —> and P — > 0. We define the diverging 
logarithmic part as our parquet variable 

A 

x = l °g-, (9) 
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It is to be noted that the higher order diagrams will con- 
tribute a factor (ux) n compared to the lowest order con- 



tribution 



Virtual processes which contribute fac- 



tors like u m x n (m ^ n) compared to the lowest order, 
are neglected in the parquet approximation^. In this 
approximation, contributions from the third order term 
in the action(equation [5] ) arc negligible compared to the 
logarithmically diverging contributions coming from the 
fourth order terms. We define the renormalized proper 
four point vertices r 4 and A4, with their bare values 



given as, 



Tl = u and = v. 



(10) 



Summing over all the leading order corrections, we get 
the following equations in terms of the parquet variable 



x, 



dT A 

d.r 
dA 4 

dx 



= --K d+l T 2 {x) 



-K d+1 T 4 (x)A 4 {x) - -K d+1 A 2 {x) (11) 



The solutions of these equations can be written as, 



r 4 = 



1 4 



1 



(12) 



It is to be noted that the parquet variable x contains 
r 4 and A 4 . Thus the set of equations [12] defines cou- 
pled equations for T 4 and A4. They need to be solved 
self-consistently to find out their dependencies on r and 
P. From r 4 and A4 thus obtained, we can calculate the 
fluctuation renormalized free energy using the relation, 



d 4 F 
dP* 
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We need to integrate (with proper boundary conditions) 
the above equation four times respect to P to get an 
expression for the free energy. Integrating the equation 
[j"3] once we get 
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Jo 

[ P'^j(T 4 + A 4 )dP' + c(r) (14) 



Here c(r) is a constant independent of P and equals to 
zero for the symmetry reason (d 3 F /dP 3 \o = 0). The 
second term in the right hand side of the above equa- 
tion takes care of the P dependence of T 4 and A 4 . To 
make our calculations simpler, we will neglect that term 
at this stage. Before doing so, we make an estimate of 
the corresponding error. The integral reads as, 
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A 4 ) 3 P' 2 dP' 



6r 



K d 



[(3r 4 + a 4 ) 

in d = 3, K d+ i(3T 4 
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Since K d+1 = ^ in d - 3, A' d+1 (3r 4 + A 4 ) - ©(lO" 2 ). 
Thus the contribution from the P-dependence of T 4 and 
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A4 to the integral (fl"5)) becomes, 

K d+lf _ . vn r P' 2 dP' 



(3r 4 + A 4 ) 3 j 
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We see that, the P-dependence of T 4 and A 4 contributes 
two order of magnitude less compared to the other terms 
in the free energy. Thus we neglect the P dependence of 
r 4 and A 4 at this stage in calculating F. Integrating two 
more times, we get, 



r)F 1 



(17) 



At this stage to have a form of F suitable to describe the 
first order transition we need to retain the P dependence 
of T 4 and A 4 . 

F = y P (rP' + i(r 4 + A 4 )P' 3 )dP' 

1 2 (r 4 (p) + A 4 (p)) 

2 4! 

- ^ P p' 4 ^(r 4 ( P ') + A 4 (p'))dP' (18) 

To evaluate the above integral, we make the following 
substitution(using equation[[S]) 
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Contribution from the integral part in the previous equa- 
tion is given by 
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Thus we have the following expression for the free energy 
at zero temperature 
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where T 4 and A 4 are P-dependent and are to be deter- 
mined from the set of equations [TT]. We notice that only 
the following combinations of T 4 and A 4 appear in all 
the calculations, 
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Here 72 is the coefficient for the quartic term in mean 
field approximation. On the other hand non-zero polar- 
ization enters into the fluctuation propagator with a cou- 
pling constant 71. In case of 72 < 0, a mean field picture 
requires an additional |P| 6 term with positive coefficient 
for the stability of the system. However a fluctuation cor- 
rected scenario can ensure stability without such a term, 
provided 71 > 0. With the above definitions, the set 
of equations [TJ becomes a single self-consistent equation 
for 71 
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(23) 



Here x contains 71 only. Solving the above equation, 72 
can be found from (eqn. |11| ) 
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Below 72 = 0, phase transition in this system will be 
first order. We are interested in the phase transition near 
72 = 0, where fluctuation effects in four point vertices are 
important. If we limit ourselves to the region |A 4 | < 3r 4 , 
the leading order behavior of 71 is same as that of T 4 and 
is given by 
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< 1 



Assuming the bare value T 4 ~ O(10), so that 
which validates perturbation theory, gives 6r!j ~ O(10 2 ). 
If we define a parameter a = jK d+ iTl, then within 
the validity regime of perturbation theory a ~ O(10 _1 ). 
Thus for an wide range of 71 e. g. 71 — (9(1) - O(10 2 ), 

For small r, 71 and 72 can be 



(20) we find |olog7i| < 



estimated as 



|2F1 

I y t 



7i 



3r° 



3r° 



(26) 



l-alog(7iP 2 ) l-2alogP 

This is essentially a non-self-consistcnt solution for 71. 
Variation of 71 at zero temperature is shown in figure 
[2]. However from the equation [24], we get, 
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where 7° is the bare value of 72 . Since the corrections due 
to self-consistency ~ log 71, the above estimate breaks 
down near P ~ exp(l/2a). Except in that regime, the 
non-self-consistent result is expected to give good result. 
Free energy at T = 0: Using the asymptotic behavior of 




FIG. 2: Asymptotic evolutions of 71 with P at T — 0. 
rameter values are 31^ = 10 and a = 0.1. 



Pa- 



Where fj, is defined as 
_ 9A' d+ i(r°) 2 



67 2 ° 



(31) 



Since >> 1, the parameter 7° > corresponds to the 
very large values of Pq. Since the scheme presented here, 
is valid for small P, we exclude this possibility in this dis- 
cussion. On the other hand, 7° < corresponds to the 
finite value of Po and there is a possibility of a first order 
transition at r = ro > 0. As r ~ P 2 , ro ~ exp - (1 + /Lt). 
If we consider the Gaussian thermal fluctuations near the 
instability point and neglect the temperature dependen- 
cies of the vertices, then the thermal corrections to r 
and P 2 ~ T 2 and Pq should vanish above a temperature 
T ~ exp y~ (1 + fx). However the parquet equations for 
the vertices and hence the form of fluctuation rcnormal- 
izcd free energy should get changed at finite temperature. 
In the next section we will discuss the finite temperature 
case in detail. 



the four point vertices [26] and using equation |21| with 

a redefinition of 7° as 7° — ^ Ti \a d+1 ■• wc can write the 
following asymptotic expression for the free energy at 
zero temperature 
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36 x 8 \1 - 2alogP, 

Here the first term ~ 0(r 2 ), hence small. The sec- 
ond term is a standard quadratic term with fluctuation 
corrections. The third term, describes the appropriate 
physics of the problem. The coefficient of the quartic 
term contains three terms. First one is a constant and 
can take either positive or negative but small values. Now 
there is a competition between the second and the third 
term. The second term tries to make the free energy 
positive while the third term tries to make it negative. 
However unless 7° is a sufficiently large negative number, 
coefficient of the quartic term is positive in the parameter 
regime of our interest. 

Possibility of a first order transition: To have a phase 
transition, the following equation must be satisfied, 

f)F 1 

■ -(r 4 + A 4 )P 3 = (29) 



rP 
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For positive r « 0, the above equation can be satisfied 
only if T4 + A4 is negative. Corresponding value of the 
P is given as, 
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IV. FLUCTUATION CORRECTIONS TO THE 
FREE ENERGY AT NON-ZERO TEMPERATURE 



FIG. 3: Diagrammatic representation of gap renormalization 
up-to one loop. 



Here we consider the system to be near the mean field 
quantum critical point with a zero temperature negative 
gap( i. e. r = — ro, ro > 0). Near ro = 0, the zero 
temperature gap attains a finite temperature value r(T) 
as follows. 
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Here the thermal fluctuations are considered up-to the 
Gaussian level, and the above expression for r(T) with- 
out any correction to four point vertices predicts a second 
order transition with a transition temperature T c ~ y'ro • 
Now we will look at the correction to the four point ver- 
tices. Using the same procedure as used for the zero 
temperature case, we now deduce the parquet equations 
for the four point vertices and hence the free energy at 
finite temperature. At T 7^ 0, fluctuation corrections to 
the four point vertices takes the form 
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in three dimension. Near the critical point, r(T) ~ r + 
Kd{u/2 + v/6)T 2 , so at finite temperature, we can define 
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our parquet variable as 



T 



In this case, fluctuation corrections to the free energy in 
terms of 71 and 72 (equation [T8]) becomes, 



(r + (u/2 + v/6)(P 2 + K d T 2 ))2 
The finite temperature parquet variable is defined as 
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In performing the above integral ^-dependence of T4 and 
A4 are neglected. For the systems near r — > limit, 
retaining only the terms lowest order in T and r, the free 
energy can be truncated as 



F = 



-rP 2 
2 

T 

- (r- 
4 L ^ 



P 4 (- 



A 4 



4! 



) 



(r + (^ + ^)(P 2 + X d T 2 )^ 



G 



(37) 



For small r and T, the third term is of the 0(T 4 ), hence 
is negligible. Thus the free energy in the leading order 
can be further truncated as 
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Following the calculations similar to that of the zero tem- 
perature case, we get the following expression for the free 
energy at the finite temperature. 
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In the above equation P = \P\ and the cubic term which 
is a result of small P expansion, does not violate the sym- 
metry of the problem. A finite temperature version of the 



equation [2fj], i.e. the asymptotic form of the self consis- 
tent equation for 71, one of the important combinations 
of the four point vertices reads, 
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Non-zero solution of the above equation tells 
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In deriving the equations at finite temperature we have 
limited ourselves in the low temperature region i.e. a 2 T 4 
is neglected compared to K d T 2 . Moreover 71 should de- 
crease at nonzero P and should go towards it's bare value 
with increasing P. Hence the part is remaining, 
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Asymptotic evolution of 71 with P for two different tem- 
perature is shown in figure [J. From equation |27| 
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FIG. 4: Asymptotic evolution of 71 with P at T ^ 0. Two 
curves are drawn at two different temperatures with V4 = 
10/3, a = 0.1 and Kd = 0.1. Upper one is at T = 1 and the 
lower one is at T = 2 in some arbitrary scale. 



Here To is the second order transition temperature for the 
mean field theory. Since near QCP, To ~ \/to, Pq is also 
~ i/tq. If we compare this results with that of the zero 
temperature case, we see that in both cases first order 
transitions is certainly weak. But in the finite tempera- 
ture case, the discontinuity in the order parameter near 
the transition point ~ ^/tq and since rg is the smallest 
scale(near mean field QCP, tq — > 0) in this system, first 
order transition as a result of the order parameter fluc- 
tuations at finite temperature near the mean field QCP 
is weaker than that of the zero temperature case. Possi- 
ble reason of such a result will be discussed in the next 
section. 



V. DISCUSSIONS 



Free energy at nonzero T: Again using the asymptotic 
behavior of the four point vertices [42] and using equation 
|39| , we can write the following asymptotic expression for 
the free energy at non zero temperature 



F 




Gap renormalization up-to one loop (figure [3]) at finite 
temperature near r = 0, tells, 



r(T) 



-r Q 



Kd 3 



6 y/P' 2 + K d T' 2 J 



T 2 . (45) 



Neglecting the term of 0(T 3 ) in the coefficient of P 2 , 
the expression for free energy becomes 

: 1 f _ , 




(46) 



In the limit P 2 » KdT 2 , it takes the form 



F 



1 K d Tl 2 

T 

24 



3aT 2 ^/Tl 



P 1 



4! 



72 V. 



(47) 



Possibility of the first order transition: The above equa- 
tion tells that, the solution of the equation dF/dP = 
will result a nonzero value of Pq ~ To, even if 7° > 0. 



Occurrence of first order transition due to coupling 
between the order parameter and soft modes is a well 
studied problem in many classical system o 12 ) 13 ) 14 ^ 
and recently in quantum phase transitions in electronic 
systems^. To capture the basic physics near a weak first 
order transitions, one needs to study the effects of fluctu- 
ations on the proper four point vertices. Bare perturba- 
tion calculations in this system show vertex corrections to 
be logarithmically singular at zero temperature in three 
spatial dimensions. We use parquet re-summation to de- 
rive recursion relations for four point vertices. Here the 
renormalized vertices crucially depends on the non-zero 
polarization. Using these relations expressions for free 
energy suitable for describing transitions at zero temper- 
ature and at finite temperature are derived. The relevant 
quantities like transition temperature and the disconti- 
nuity in the order parameter at the transition point turn 
out to be small but finite. We mainly concentrate on the 
phenomena near a quantum critical point predicted by a 
mean field theory. We found stronger possibility of first 
order transitions at T = than at any finite tempera- 
ture transition near QCP which is in accord with recent 
experimental finding. It is found that a finite polariza- 
tion tries to suppress fluctuations. In such a scenario, we 
found that exploiting the logarithmic singularity in the 
vertex corrections, the zero temperature system is able 
to gain a finite possibility to undergo a first order tran- 
sition near mean field QCP while the finite temperature 
one fails to do so. Moreover unlike the standard RG ap- 
proach, which considers first order transition as just the 
inability of the system to reach an unstable fixed point, 
the present approach makes qualitative predictions about 
the magnitude of the discontinuity in order-parameter 
near the transition point. Similar phenomena occurs in 
low T c itinerant magnets because of the coupling between 
order parameter and other soft modest. The effects 
of fluctuations in the discontinuity in order-parameter 
has been estimated earlier using loop corrections to the 
free energy, hence an non-analytic Landau expansion. In 
those works the occurrence of first order transition at 
low temperature is concluded to be generic. We look at 
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the similar problem in a different physical system and 
emphasize the importance of competition between fluc- 
tuations and nonzero polarization and the crucial role 
played by marginal dimensionality. We have considered 
the constant pressure case only where strain effects are 
completely integrated out. One can also consider a sit- 
uation where system volume is constant and a proper 
quantum generalization of such case^ should certainly 
be addressed. 
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